Submitted to ApJ Letters 

Preprint typeset using I£T£jX style emulateapj v. 14/09/00 



EVIDENCE FOR SOLAR-LIKE OSCILLATIONS IN BETA HYDRI 
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ABSTRACT 

We have made a clear detection of excess power, providing strong evidence for solar-like oscillations in the 
G2 subgiant (3 Hyi. We observed this star over five nights with the UCLES echelle spectrograph on the 3.9-m 
Anglo-Australian Telescope, using an iodine absorption cell as a velocity reference. The time series of 1196 
velocity measurements shows an rms scatter of 3.30ms" 1 , and the mean noise level in the amplitude spectrum at 
frequencies above 0.5 mHz is 0. 1 1 m s" 1 . We see a clear excess of power centred at 1 .0 mHz, with peak amplitudes 
of about 0.5 m s" 1 , in agreement with expectations for this star. Fitting the asymptotic relation to the power 
spectrum indicates the most likely value for the large separation is 56.2 /iHz, also in good agreement with the 
known properties of (3 Hyi. 

Subject headings: stars: individual (J3 Hyi) — stars: oscillations — techniques: radial velocities 



1. INTRODUCTION 

The search for solar-like oscillations in other stars has been 
long and difficult. Observers have mostly concentrated on three 
stars: Procyon (a CMi), r\ Boo and a Cen A. Reviews of those 
efforts have been given by Brown & Gilliland (1994), Kjeld- 
sen & Bedding (1995, hereafter KB95), Heasley et al. (1996) 
and Bedding & Kjeldsen (1998). More recently, Kjeldsen et al. 
(1999) measured B aimer- line equivalent widths in a Cen A and 
set an upper limit on oscillation amplitudes of only 1 .4 times 
solar, with tentative evidence for p-mode structure. More re- 
cently still, measurements of velocity variations in Procyon by 
Martic et al. (1999, see also Barban et al. 1999), showed very 
good evidence for oscillations with peak amplitudes of about 
0.5 ms" 1 and frequencies centred at about 1 mHz. 

Here, we report the clear detection of excess power, provid- 
ing evidence for oscillations in the G2 subgiant (3 Hydri (HR 98, 
V = 2.80, G2 IV). This star is the closest G-type subgiant, with 
luminosity 3.5 L©, mass 1.1M© and age about 6.7 Gy (Dravins, 
Lindegren, & VandenBerg 1998). It has received less attention 
than the three stars mentioned above, presumably due to its ex- 
treme southerly declination (-77°). An attempt to measure os- 
cillations in (3 Hyi in radial velocity was made by Edmonds & 
Cram (1995) and gave upper limits on the strongest modes of 
1.5 to 2.0 ms" 1 , consistent with our detection. 

2. OBSERVATIONS AND DATA REDUCTION 



55±0.1°C and imprints a rich forest of molecular iodine ab- 
sorption lines from 500 nm to 600 nm directly on the inci- 
dent starlight. Echelle spectra were recorded with the MITLL 
2kx4k 15/im pixel CCD, denoted MITLL2a, which covered 
the wavelength range 470-880 nm. 

Exposure times were typically 60 s, with a dead-time of 55 s 
between exposures (using "FAST" readout). The signal-to- 
noise ratio for most spectra was in the range 200 to 400, de- 
pending on the seeing and extinction. In total, 1196 spectra 
were collected, with the following distribution over the five 
nights: 59, 166, 301, 325 and 345. The first two nights were 
affected by poor weather and technical problems. 

Extraction of radial velocities from the echelle spectra fol- 
lowed the method described by Butler et al. (1996). This in- 
volved using the embedded iodine lines both as a wavelength 
reference and also to recover the spectrograph point-spread- 
function. Essential to this process were template spectra taken 
of (3 Hyi with the iodine cell removed from the beam, and of 
the iodine cell itself superimposed on a rapidly rotating B-type 
star. 

The resulting velocity measurements for f3 Hyi are shown in 
the lower panel of Fig. 1. They have been corrected to the 
solar system barycentre, as described by Butler et al. (1996). 
No other corrections, decorrelation or high-pass filtering have 
been applied: the measurements are exactly as they emerged 
from the pipeline processing. Note that the velocity measure- 
ments are relative to the velocity of the star when the template 
was taken. The slow variations in velocity within and between 
nights are due to uncorrected instrumental drifts. The rms scat- 
ter of these measurements is 3.30 m s" 1 . 

Uncertainties for the velocity measurements were estimated 
from residuals in the fitting procedure and are shown in the up- 
per panel of Fig. 1. Most lie in the range 2.5-4 ms" 1 , and fall 
gradually during each night as the target rises in the sky. 



The observations were made over five nights (2000 June 1 1- 
15). We used the University College London Echelle Spec- 
trograph (UCLES) at the coude focus of the 3.9-m Anglo- 
Australian Telescope (AAT) at Siding Spring Observatory, Aus- 
tralia. To produce high-precision velocity measurements, the 
star was observed through an iodine absorption cell mounted 
directly in the telescope beam, immediately behind the spec- 
trograph entrance slit. The cell is temperature-stabilized at 
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FIG. 1 . — Velocity measurements of /3 Hyi obtained with the AAT (lower panel) and the corresponding uncertainties (upper panel). 



Observations of (3 Hyi were also made using the CORALIE 
spectrograph on the 1.2-m Leonard Euler Swiss telescope at La 
Silla Observatory in Chile (Queloz et al. 2000). The precision 
of those measurements was poorer than those from UCLES, 
presumably due - at least in part - to the smaller telescope aper- 
ture. The CORALIE data are not included in this Letter; their 
analysis is postponed to a future paper. 

3. TIME SERIES ANALYSIS 

The amplitude spectrum of the velocity time series was cal- 
culated as a weighted least-squares fit of sinusoids (Frandsen 
et al. 1995; Arentoft et al. 1998), with a weight being assigned 
to each point according its uncertainty estimate. Figure 2 shows 
the resulting power spectrum. There is a striking excess of 
power around 1 mHz which is the clear signature of solar-like 
oscillations. As discussed below, the frequency and amplitude 
of this excess power are in excellent agreement with expecta- 
tions. We also note that the excess is apparent in the power 
spectra of individual nights. 

Typically for such a power spectrum, the noise has two com- 
ponents: 

1. At high frequencies it is flat (i.e., white), indicative of 
the Poisson statistics of photon noise. The mean noise 
level in the amplitude spectrum in the range 3-5 mHz is 
0.11 ms" 1 . Since this is based on 1196 measurements, 
we can calculate (e.g., KB95) that the velocity precision 
on the corresponding timescales is 2.2 ms" 1 . 

2. Towards the lowest frequencies, we see rising power that 
arises from the slow drifts mentioned above. A logarith- 
mic plot shows that, as expected for instrumental drift, 
this noise goes inversely with frequency in the ampli- 
tude spectrum (and inversely with frequency squared in 
power). It has a value at 0.1 mHz of 0.3 m s" 1 . 

The white line in Fig. 2 shows the combined noise level, and 
we see that the contribution from 1// noise is neglible above 
about 0.5 mHz. 

4. DISCUSSION 

4.1. Oscillation frequencies 

Mode frequencies for low-degree oscillations in the Sun are 
reasonably well approximated by the asymptotic relation: 

v(n,l) = Av(n + \l + e)-l(l+\)D a . (1) 
Here n and / are integers which define the radial order and an- 
gular degree of the mode, respectively; Av (the so-called large 
separation) reflects the average stellar density, Do is sensitive 
to the sound speed near the core and s is sensitive to the sur- 
face layers. Note that D is Sv /6, where Sv Q is the so-called 
small frequency separation between adjacent modes with / = 
and / = 2. 



A similar relation to equation (1) is expected for other 
solar-like stars, although there may be significant deviations in 
more evolved stars. Models of the G subgiant 77 Boo suggest 
that modes with 1=1 undergo 'avoided crossings,' in which 
their frequencies are shifted from their usual regular spacing 
by effects of gravity modes in the stellar core (Christensen- 
Dalsgaard, Bedding, & Kjeldsen 1995). Some evidence for 
this effect was seen in the proposed detection of oscillations 
in rj Boo by Kjeldsen et al. (1995). 

In attempting to find peaks in our power spectrum match- 
ing the asymptotic relation, we were severely hampered by the 
single-site window function. As is well known, daily gaps in a 
time series produce aliases in the power spectrum at spacings 
±11.57 /iHz which are difficult to disentangle from the genuine 
peaks. Various methods have been discussed in the literature for 
searching for a regular series of peaks, such as autocorrelation, 
comb response and histograms of frequencies (e.g., Gilliland 
et al. 1993; Kjeldsen et al. 1995; Mosser et al. 1998; Barban 
et al. 1999; Martic et al. 1999). Here we use a type of comb 
analysis, in which we calculated a response function for all sen- 
sible values of Av, Do and e. 

Since we are searching for mode structure in the region of 
excess power, it is convenient to rewrite equation (1) as 

v{ri ,/) = vq + Av(ri + ±/)-/(/+l)D . (2) 

Here, vq is the frequency of a radial (/ = 0) mode in the region 
of maximum power, for which n' = 0. Thus, v replaces e as the 
parameter for the absolute position of the comb. 

The comb response function was obtained as follows. We 
first thresholded the power spectrum at the noise level. In other 
words, any points less than the noise level were set at this level. 
For each triplet {Av, Do, vq) we then measured the (modified) 
power spectrum at each of the frequencies predicted by equa- 
tion (2). We summed these numbers to produce a response 
value which is a measure of the goodness-of-fit of that triplet. 
We limited the sum to / = 0,1,2, and n' = -4,-3, ... ,4, with vq 
centred at 1000 ^Hz. The restriction on n' was made for two 
reasons: as in the Sun, we only see modes excited to observ- 
able amplitudes over a restricted range of frequencies and, also 
as in the Sun, we expect deviations from equation (2) over large 
ranges in frequency. 

From this analysis we identified the parameters that max- 
imised the comb response, searching over the ranges Av = 45- 
80 /iHz and Do = 0.5-2.0 /iHz. These ranges encompass the 
values that would be expected for this star. For vo, we exam- 
ined a range of width Av centred on 1000 ^Hz. The follow- 
ing two triplets of (A^, Do, vq) gave the best responses: (56.2, 
0.83, 1030.1) ^Hz (solution A) and (60.3, 0.75, 1005.3) ^Hz 
(solution B). The frequencies corresponding to these solutions 
are shown in Fig. 3, overlaid on the measured power spectrum. 
The distribution as a function of frequency of matches between 
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FIG. 2. — Power spectrum of the AAT velocity measurement of fi Hyi. The white line shows a two-component noise model (see Sec. 3). 



calculated and observed peaks is not the same for the two so- 
lutions. The matches in Solution B occur mostly at the lower 
frequencies, while those for Solution A cover the full range of 
excess power. This suggests that Solution A is more likely to 
be the correct one. 

How do these two solutions agree with expectations for 
(3 Hyi? The large separation Av of a star scales approximately 
as the square root of density. Extrapolating from the solar case 
(Av = 135.0/iHz) using parameters for (3 Hyi of L = 3.5 L©, 
T eS = 5800 K and M = 1 . 1 M (see Dravins et al. 1998, and ref- 
erences therein) gives Av = 56.0 yuHz. Given that we searched 
over the entire range 45-80 /iHz, the excellent agreement be- 
tween observation and theory is very encouraging, especially 
for solution A. 

There is no simple way to estimate Do for a star as evolved as 
(3 Hyi, since this parameter is very sensitive to evolution. Ex- 
trapolating the grid calculated by Christensen-Dalsgaard (1993) 
leads us to expect Do in the range 0.5 to 1.0, which is consis- 
tent with both solutions. The value of vq (or, equivalently, e) 
depends on details of the surface layers, which also requires 
construction of models specifically for (3 Hyi. 

We note that simulations in which input frequencies obeyed 
the asymptotic relation precisely showed a more definite (and 
unambiguous) peak in the comb response than did the real data. 
We conclude that the oscillation frequencies in (3 Hyi have sig- 
nificant departures from the asymptotic relation. This is of great 
astrophysical interest and not unexpected, but it does make it 
very difficult to extract the correct frequencies from our single- 
site data. 

At the suggestion of the referee, we examined the comb re- 
sponse of random power spectra that were generated by mul- 
tiplying white noise by an envelope similar to the observed 
power excess. For these spectra, the comb response sometimes 
showed peaks similar to those from the real data, at various val- 
ues of Av. The comb analysis therefore does not prove that 
the excess power is due to individual modes. However, given 
(i) the agreement between observations and expectations in all 
respects - position and amplitude of excess power (see below), 
and the best-fit value for Av - and (ii) the absence of any the- 
oretical or observational reason to attribute the power excess 
to other sources, we consider that solar-like oscillations are the 
most likely explanation for the observed power excess. 

4.2. Oscillation amplitudes 

The strongest peaks in the amplitude spectrum of (3 Hyi 
(square root of power) reach about 0.6 m s" 1 . However, these 
are likely to have been strengthened significantly by construc- 
tive interference with noise peaks. As stressed by KB95 (Ap- 
pendix A. 2), the effects of the noise must be taken into account 
when estimating the amplitude of the underlying signal. To do 
this, we have generated simulated time series consisting of arti- 
ficial signal plus noise. We conclude that the underlying oscil- 
lations have peaks of about 0.5 m s" 1 . 

Solar-like oscillations are excited by convection and the ex- 



pected amplitudes have been estimated using theoretical mod- 
els. Based on models by Christensen-Dalsgaard & Frandsen 
(1983), KB95 suggested that amplitudes in velocity should 
scale as L/M. More recent calculations by Houdek et al. (1999) 
confirm this scaling relation, at least for stars with near-solar ef- 
fective temperatures. For (3 Hyi, the implied amplitude is about 
3.2 times solar, which is 0.7 to 0.8 m s" 1 (based on the strongest 
few peaks - see KB95). The observed amplitudes are therefore 
consistent with expectations. We also note that the frequency of 
excess power (1 mHz) is in excellent agreement with the value 
expected from scaling the acoustic cutoff from the solar case 
(Brown et al. 1991, KB95). 

5. CONCLUSION 

Our observations of f3 Hyi show an obvious excess of power, 
clearly separated from the 1/f noise, and with a position and 
amplitude that are in agreement with expectations. Although 
hampered by the single-site window, a comb analysis shows 
evidence for approximate regularity in the peaks at the spac- 
ing expected from asymptotic theory. There seem to be signifi- 
cant departures from regularity, perhaps indicating mode shifts 
from avoided crossings. We hope that further analysis along 
one or more of the following lines will allow us to explore fur- 
ther the oscillation spectrum of (3 Hyi: including the CORALIE 
observations which, despite their lower precision, could use- 
fully improve the spectral window; reducing the noise of the 
velocity measurements (UCLES and CORALIE) by decorre- 
lating against external parameters (e.g., Gilliland et al. 1991); 
combining these results with an analysis of equivalent-width 
variations of strong lines in the spectral region uncontaminated 
by iodine lines (e.g., Ha); and using theoretical models to cal- 
culate expected shifts due to avoided crossings, to help identify 
the affected modes. 

The clear strength of the technique used in these observa- 
tions is its ability to precisely calibrate spectrograph variations 
at the timescales of primary interest. Although the best ra- 
dial velocity precisions achieved in the long term from this 
technique are 3-4 m s" 1 , our data clearly demonstrate that the 
precision at frequencies around 1 mHz is significantly better 
(2.2 ms -1 ). Our results provide valuable confirmation that os- 
cillations in solar-like stars really do have the amplitudes that 
we have been led to expect by extrapolating from the Sun. This 
bodes extremely well for success of space missions such as 
MOST (Matthews et al. 2000), MONS (Kjeldsen, Bedding, & 
Christensen-Dalsgaard 2000) and COROT (Baglin et al. 1998), 
which will provide photometric data of high quality for a wide 
sample of stars. 
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(RPB), SUN Microsystems, and the Danish Natural Science 
Research Council and the Danish National Research Founda- 
tion through its establishment of the Theoretical Astrophysics 
Center (HK). 
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FIG. 3. — Close-up of the power spectrum of /3 Hyi, with dashed lines showing the frequencies given by equation (2) for solutions A (upper panel) and B (lower 
panel). The dotted lines show the ±l/day aliases of those frequencies. 
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